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ABSTIBACr 


Filter aircraft often encounter handling difficulties during 
the crucial apprcach phase of the landing maneuver. Ihis manifests 
itself as a significant dip in tlie flight path angle before 
settling to the desired command levels. Ihis abnormal trend can 
generate considerable confusion in the pilot during this vital 
flight phase. At present, tlie pilots attempt to cope with this 
airplane behavior by constant Juggling between the throttle and 
the stick. It may be pointed out tliat this pilot maneuver is 
rather damandingi especially in view of tlie large lags in the 
response of the mgine. 

In the {sresent investigation an attempt has been made to 
understand the root cause of the problem through a parametric 
study. Ihe study suggests that the presence, in the right lialf 
s-plane, of the zero (~1/T ) in the flight path angle - to - 

hi 

elevator transfer funoticxi, brings about the undesireable flight 
path response during landing approach. 

Next an attenipt has been made to synthesise a MIMO 
control system utilising the airbrakes, in order to ensure that 
this zero stays in the left half s - plane. A classical approach 
to MIMO systems is adapted to the design problem at hand. This 
approach utilises the closure of llie feedback loops in a 
seqpentlal manner. It also allows for the placement of the poles 
as well as the zeros of the transfer function of interest by a 
simple root locus design procedure. 



The results clearly demonstrate tlie efficacy of the proposed 
autopilot in overcoming the problem of the flight path dip during 
the landing phase as well as statically decoupling the airspeed 
from the flight path. 
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CHAPTER 1 
INl'RUUUCTI ON 


' The landing approach represents a ratiter iaportnnt phase 
in an airplane flight mission. Tlie suocer.sful ■.■(«¥> lotion of the 
mission requires much greater precision in airplane handling 
daring this crucial phase. The inportance of l-Jiis becomes evident 
from the large number of failures and associated loss of life that 
nay be attributed to inadequate flying qua li hies during this 
terminal phase. 

The acceptable handling quaJities for the approach phase 

f 

demand control over two flight parameters - hlie frtrwani air.apeed 
as well as the rate of descent. Tiie pilot's twindling hank involved 

here is particularly moire demanding heoaiine of tlie coupling 

« 

invariably present between tlie vertical sink rate atd the 
airspeed. Evidently, this significantly adds to tlie pilot's 
workload. Heedless to say, in the present era of modem hi|^ 
performance airplanes there is little scope for pilot error even 

m 

in the face of such excessive demands, more so in this terminal 
phase. 

For the twin control of tlin sink ralx* and tlie airspeed, only 
a single aerodynamic surface is available - tl>> elevator. The use 
of the throttle enables tlie niri >19110 designer to provide for tlie 
second control normally neederl. To afhievi- proper ncordination 
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between these two flight peranieLer.s tlie pilot, mnhen use of tlio 
throttle on ni 8£:-atid-wiM»i-rieeii>«l I'Mf:)::. Kf»pl.iitifi* uly this 
furiotion does not normnlly involve- ••ru-fi-Undn t ■•riiili>i! 

eotuntiori. 


In genernl, it is eas'.y r'ii:»ijv uiifdi-ftoL'ii I IvMKlJiug 

in iion-tjsrmijiel pimses; li*.«Hevfr, 'Iiirin?! i.t-** luxuliiig HitTori'Hi. 
there bctiiig ewtoseiivA Mrp* l.-ir.; Jn Ui* 

difficulties are encounLi're.1. (hiA i-*;vi cmly 1ni.<»giiio tiie probleiis 
while landing upon an airornft cnrtior nt sea. c-n|>e'oiaiJy duidjig a 
stom. Hence providing tlie pilot with a suitrilde control mecliaiiism 
without such a time lag proLilem SEipears tu lie a mjst. Uveri duiing 
the less critical phases it would lie worthwhile bo have dual 
siiultaneous control over the rate of descent and tlie airsprjed tu 
reduce the pilot's overall workload aivl tlus indirectly enliaiioe 
nission success rates. 

The present controller designs ll»P] i-ciitre aruuiid tliu idea 
of using feedback actuation of tlie thi-ottle ‘and tlie elevator 
simltaneoualy. Tliese ss^tems suffei' fixm tie throttle tiuo-lag 
problem. Besides, the thrust needs to be adjusted to re-establiuh 
equilibrium each tine there is a disturbance or when the pilot 
wants to change either of the two fliglit paramoiers. This causes 
the wear and teat of the power plant to inorcuse significantly 
121 . Thus, for better liandliiig ami i-edui'ixl riigiiitt wear, .an 
alternative actuaticn system is rr^qult-ed in tlie iongitudliia] 


chaiuiel. ' 
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The airbrake is selccLul lieri! tci ptuvldv rut- tlie tvJdLbioiiel 
control surface. ]b nay seem uii uiiuuuveaitiMiiul ^'liovre, Liut it 
offers several advantages over tlie tlii-uttlt; rr’;i. FirntJy, tlie 
control being exercised through an aex-odyriamio sui'fuue, it ItuidJy 
suffers from tlie time-lag profcilen asctociated with tlie thiottie. 
Secondly, it reduces tlie wear and tem* nn the engiiu', imw 
operating at a constant tlirnttln netting. [■' inaJJy. siiiCH* l.lie 
airbrake is d^loyed, the tlirottJe seti.ing is gnlng to be higiier 
than it would be during a conventional hindirig. Tlut turns out to 
be an added advantage for a suddon abortion of landing — if 
forced — as the pilot can siifly retract tlie airbrake to get an 
alnost instantaneous boost in the airspeed. Tlie gu-around is 
further facilitated through faster engine response by virtue of 
the higher throttle setting now Ui operation. JL is felt Uiab 
these advantages may make Uie airbralio a particularly 
attractive choice for the design of tlie flight, pal.h coiiLruiler for 
the approach phase of tin ixidem high performance fighter 
airplanes. 
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CIIAI’IRK rt 
I‘R(MI1.FM ANAl.ySrS 


For a better assesssieiit of ihu probl^.u oL' ivsgiiJatiiig the 
flight path during appi-oach, it in d'-sii-riidf^ |j> ik'Vfslcfi n i.-lr'ai' 
UTKlerstanding of tlie various Lerun iiivulvi-fd in tlie governing 
transfer function of tlie airplane. This wi'niid enable us to 
identify the sore isportaiit systeai parameters and hence tlie 
associated stability derivatives. 


We therefore carefully examine a Lypicral transfer ftviction 
for a fighter aircraft during tlv^ tj^tsiinal pliase umjer 
consideration 131: 


y(s) 

«,<»> 

where 


to' (s + 4 + 1/1',^.,) 


.Cs*+ zr. *»> o 4 i.>* JLn*i <7 

^ •!* I'll I'll |M| 


r - flight path angle (rad; 

6 = elevator defleutiui (rail) 

0 

s = Laplacian varisjibld 
ky = constant gain factor 

= short period danping ratio 
CO = short period natural freguenoy (rad/s) 

®P 

C . = plugoid daigiing faoLor 

pn 

CO . = pliugoid natural fi-ejiioiicy (rad/s) 

pn 

(J/T^^), Cl/r,,,,) and (1/T,^j^; artr numurutor LIbic ouiiuLuntc with 
1) T„ > 0 ..Kl • M 

11 ) I T^. I I T,^ I nr I T,„ I 
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It is evident that the above transfer funubioii lies : 
i) zeros at 

U) voles at Ct-t.p i J 

t«ph * > J-'’ph J 

The iiportance of these s-plane Juontjunu cannot be over 
esphasised since the airplane flight path response during any 
flight phase is dictated by their iiistjuitmieous vaLuos. As is well 
knoMfi, the above parameters tindergo significant variations as the 
airplane goes through changes of flight regime and maneuvers. 

It is found that increasing t)» Mach number causes the two 
extreme zeros ( and to drift further a^t while 

the third. (-l/T^^). lying between these two and rather close to 
the imaginary axis, moves to ttie left. Inversely, during the 
oritioal phase of approach, when Hach numbers are low. this middle 
zero shifts to the right* getting invariably puslied even into the 
right-half plane. This appears to be at tlie root of the control 
problw. This fact can be better illustrated through the example 
given below. 

b 

Exaiqple : Typical data for a fighter ocnfiguration : 

■MM— = 9000 kg length = 12 m 

? wing area = 40 sq m m.a.c. = 7 m 

' flight oondition : Sea level. M = 0.2. r = -2 d^ 



G 


rCs) 

Vs) 

0 . WKKs - M .CW.S2> . f 5S ^ A . ) . < s - 29W > 

[s*+ 7.(0.88X1.36)3 h ||r.*i ;!(H nxy.JfJ^.n I (Ci.lH)*') 

The fllfSlit r>al.h reBpunrwf l.ri r elr^veliit- iiit.<iit. jit rIkiuii in 

Fjg.l. The Pjrst thine himl rs-in l'<* iinLli-'fl ff'ei l.lt^ is 

tliat it exhlhlt.3 a non-minimnm plmne htf.*h»ivifiiir i.n.. i.i>“ flighl. 

path angle first dips to one side and then settles at a value of 

the opposite sign. This, in i Itself , is ^aI^Pio1ettt cause 

serious problens in airplane handling since the pilot feels "out 

of phase" sith the airidane. The second in(H;irtant thing to be 

noted is that a positive elevator in|.njt onuses a positive steady 

state value of flight path. This has rhvimir. inidesirahle flying 

isplications as now pushing forwnrrl rm Uic piiot'n sl.iok itnutHS 

the airplane to settle in a si-eedy olisfci - in direct Cf.<itr8st to 

» 

the mind'-set of tlie pilot for tlie (il.her Plight pluses at liiglier 
Mach numbers. Evidently this ean leaJ f.o sliustictis-i op confusion 
■and occasionly even disaster. For an understanding of this 
observed behaviour diaring the approsicii phast;, we consider tho 
transfer ftinction in a more convenient fom : 
y(8) (s a).(s + b>.(s i- c> 

^^(s> [»*+ Ps + q]-[K*+ Is + ml 

where : 

b = VT^ 

o = 




( Bap )aujULJ ag 


OOZ- 


8 


p 

q 


2.C .o 

•p •p 


Cl> 


ap 


i = 2.^ ^ 

ph ph 



Hote that ki' is taken equal to 1 for siat'licity, (lowever It 
Ims no inplications as regards tlv* validity of tlie analysis 
whatsoever. 


It nay be repeated that the overall flight path response 
depends upon the above paranater values as well as initial 
conditicns. However« for the purpose of analysis, we assume 
initial conditions to be zero . This would enable us to ' highlight 
the intrinsic behaviour under the influence of the iifiortant 
pajrameters. Again, focus ing upon the same objective, we take the 
poles, which characterise the airplane short period and phugoid 
modes, to be fixed while the zeros are varied. 

There being three variable parameters a,b and c, these are 
varied one at a time while the other two ax-e Itopt constant, 
furthermore, since experience shows that tlie two extreme zeros 
(-1/T^ and -1/T^^) move apart in unison, the number of parameters 
ge.t8 effectively reduced to only two — a and b (or c>. 

It is interesting to point out blet to examine the changes in 
the nature of the response, Jt Is licit, iiocxifssary t.o solve bhn 
governing differential equation for the flight path. In fact it 
would suffice to break the transfer fuiiutloii- into partial 
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fractions in ihs form 
rCa^ A.Cs+B:> CCs+OD 

S 4. 

6^Cs3 [s*+ ps + q] Cs*+ Is + ml 

and simply focus upon ihs nature of the constants A»B»C and D and 
study how they change when the parameters a»b and c are varied. 
Needless to say» relatively large changes In the values of these 
constants would Indicate sudden changes In the nature of response. 
The same could be concluded from sudden changes In the signs of 
the emstants. Of the t%«o modes* the phugold being more InqxMrtant 
for the flight path* the associated constants C and D require 
closer scrutiny. 

Results of the Parameter - Study t 

The study undertaken here shows that of the three zeros* 

the one at CHL/T. 3 Is the most significant one as the changes In 

hi 

Its value can bring about radical changes In the nature of the 
longitudinal response. Of the two modes* however* It Is the 
phugold vriilch suffers more drastic behavioural changes while the 
short period one remains virtually unaffected. 

Table 1 shows the results of the parameter - study. In the 

first part* the constant Cl/T.^^ Is varied from -*0.0 to 40.5* the 

hi 

total variation being 12509( around Its nominal value* mAille the 
other two parameters and are held constant at 

their nominal values of -8. 3 and +4. S respectively. 



TAWL1-: 1 . 

RESULTS OF THE FARAHI-irER S'lUliV 


Y / ^ Nnme^ohor 

TiIkii h IVi iiirl 



Hiniiir:iU'>t' 

tlinift'iil.i.ii' 

(9fa)(sH-b)(sfc) 

ACs + B) 

C(s + D) 

(s44.2)(8-2.3)(s-«.G) 

L0.2l.Ks \ 2.LVI) 

•■3.28(8 - 0.2ii) 

Csf4.2}(8-2.3)(s-0.2) 

7.01(8 4 2.17) 

■■6.91(« - H.12> 

C9f4.2}(s-2.3)Cs-0.1> 

7.12Ln \ 2.;^) 

•R.12(b - U.05) 

(sl-4.2)(s-2.3}(s-0.08) 

6.?lR<s + 

-.'l.«6(s - 0.03) 

C»f4.2)C8-2.3)Cs-0.05> 

R.72^s + 2.27) 

•5.72(8 - 0.01) 

(8f4.2)Cs-2.3)(s) 

e.•n(^• + 2.:n > 

-5.33(8 + 0.04) 

(84-4 . 2) C s-2 . 3>(8f0.05) 

ri.9'K« + 

-4.93(S + 0.00) 

<8f4.2)Cs-2.3)(&t0.08) 

5.69(8 1 2.39) 

-4.69(8 4 0.J3) 

(8f4.2)(8-2.3)(8+0.1) 

ft.r.4(n -4 2.'12) 

-4..‘VWs 4 0.18) 

(8+4 . 2) ( 8-2 . 3 )(s+0 . 2) 

4.74(8 »■ 2.G!i) 

3.74(8 4 0.32> 

(sf-4.2>(s-2.3)Csf0.5) 

2.37(8 •» 3.M) 

...... 
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To sUidy tlie effect, of VHX'ying Uio ccniftl-nnl. ) npun tlif* 

nature of tlie two Ic^igiturlinal nwlr*:: w l.ho fxtfiiii. of 

variat.ione in tlie values nf Uir* ;:ei-n:: of tli- i.w<« :'li<^iiil. Lli'-ii- 

nominal values. It Is found ttmt lli«# ''.i.a-o of Hie ytK^rt pevl»jd node 

varies only 65% around nominal while tho plnujoid xret-n varies '/1V5% 
around nominal. Furtlieniore, the sif^i of (l/’f,,, ) of epocial 

significance here an it dictates liic very nature of tlvi plnignid 

mode. For ttie case O/Tj^^ ) W. lire |.•llln*^lill ''xhiliilj*. n 

non-fflinimum pliase behaviour which is chai‘Etf'teii::e«.i by a sorn in 

the riglit-lialf plane. On tlie otlicr liaiid. wlieii this 

abnormality disappears. 

Nexij we focus on studying the influonce nf varying tlm trims 
<VTj^^) and (VT^g). It moy lie added tliat the verial.ions of Llicne 
two zeros lieppen to have symmutty arouml tlie imaginary axin an 

tliey move apart. Tlie term > Is »r.iw Imld conulaiit al Itn 

nominal value of i^.08. 

In this case, partial ft action dencat>uoitiori allows tliat 
except for a variation in the cuRiotant gain, tlx# 'nature' of 
neither mode changes significantly for tlie zeros rcnmui virtually 
unchanged. As a result, merely tlieir steady state values undergo 
variations. 

In tlie final phase, IJv< effci't nf vniyini! I.I 10 cutii:loiil.R 
(l/Tj^) andd/I'j^j,) is exaiilncil wlu'ii llw niiui “f t l/f, ^ ) in 
rcv«rrw<I, liil.rrr!:l.hii;ly, l.ln* "f lli** I ihrni:i| 



TABLE 1. CroiitU. J 
RESULTS (}F niE l*ARAMETEk SI in>Y 


r / 6m Mamfrabor 

SlKiit [V'l-iod 

riaignid 


Nuaierator 

Munerator 

(sfa}C»^})(84«} 

A(8 + B) 

C(s + D) 

CB-0.08)(Bf4.2)(8-2.3) 

8.S6(s + 2.25> 

-5.96(8 - 0.03) 

(8-0 .08 } (8+5 .0)(8-3 .0} 

10.20(8 + 2.28) 

-9.20(8 - 0.04) 

(8-« .06 ) (8+8 .0)(8-6 .0) 

30.22(8 •+ 2.20) 

-28.22(8 - 0.04) 

(8-0 .06) (8+10 .0}(8-8.0} 

55.65(8 +■ 2.28) 

-54.65(8 - 0.04) 

(8-0.08)(8+12.0)(8-12.0) 

88.36(8 + 2.31) 

-67.36(8 - 0.04) 

(840.08)(8M.2)(8-2.3) 

.5.96(8 -+ 2.38) 

-4.69(8 + 0.13) 

(8+0.08 )(fl^5 .0X8-3 .0) 

8.32(8 •+ 2.43) 

-7.32(8 +■ 0.13) 

(8+0.08}(8+8.0)(8-6.0) 

24.59(8 ■+ 2.45) 

-23.59(8 + 0.12)- 

(8+0.08)(8+]0.0}(8-8.0) 

45.33(8 •+ 2.46) 

-44.33(2 4 0.12) 

(840.06)(8+12.0}(8-. 12} 

71.88(8 •+ 2.46) 

-70.99(8 4- 0.12) 







observed esrller remain unchanc^ 

■ Thus it may be ooncluded that it is the si|;n as well as the 
magnitude of <1/T > that really decides the nature of the phugold 

hi 

mode. This mode beine more important to the flight path response^ 
we arrive at a general and rather far-reaching observation that 
the location of the zero C'-l/T|^^> is crucial in dictating the 
nature of the airplane's flight path behaviour. 

Before going on to the design of the controller, it is 
worthwhile to try and understand the physics behind the strange 
behaviour of the flight path response. It becomes quite clear when 
one realises that flight path response is largely dependent upon 
the interplay of changes in lift and drag when the elevator Is 
deflected. In this context it may be pointed out that during the 
approach phase, the airplane is operating on the left side of the 
'drag vs. speed' curve . 

Now we consider . what happens when the elevator is deflected 
downward. The resulting nose-down moment causes decrease In the 
aiq;le of attack with the attendant reduction in drag. This causes 
the airplane to accelerate as the airplane is operating at a speed 
below the minimum - drag speed. The resulting increase in speed 
leacta to further reduction in drag which in turn causes the 
airplane to accelerate even more. The large gain In speed mads 
possible through elevator deflection is accompanied by ' a 
corresponding hike in the lift, causing the flight path angle to 
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shift. to a positive value. Thh: l.iimJ is lilglily 

undesirable, it boing tin smiii-i* nf coiifii:«i>.4i lii l.ln iiilol. ;iir. 
pointed out earlier. It in l.hnr tlml in onier tn (■•mliol 

the flight path at loc? Mach nuailwsrs:, it. Is ini^eretivt* tu 
incorporate an (AklitioiiaJ »^>ri«l fmit.tril. TIk' ■ hi i-f.-iiL I'l.sr-tjc.e tUis 
involves providing such ci4itiiil |•ll.ill'^l‘ mnini.il ly op t.linmgli nn 
autopilot. 


Ho now direct our atientiim t>i tla* "nil inpotlant" sevj given 
by the cotistaiiL (VT^^). This nay be expresstxl as IdJ : 





Zu '' “cf.e-2u> 




26e-»<v / f'6e-Z«> 


TIus this zero can be varied by cliunging tlie stuiiility 
derivatives and X^. Here this speed I'endtiaok oontroJ oilmlively 
augments tin derivative . Tills in tnin rkingfs tlie nngnil.uilci uu 
well tin sign of the zero (1/1',^^). This is precisely what is 
atteupied in tin design of tin fligliL patli c*r«itrolier. 



CHAPTER 3 


CUNIRCILLER SYNl'Ill^SLS 









17 


wtere : 


A = 


a a s 

<1 ta 19 


*1 


■S* «2. 


a_ 


"91 ®92 ®9b 

Fron (l; and (2) we get 

A.X = B.(O^.X^ - G.X) 
c c 

=> (A + BG).X = B.G .X 

c c 


; B .. 


f f 


21 


91 


19 . 




9 . 9 . 


9X 


( 


The correEponding characteristic {.•olynonjal utf tlic closed loop 
system is 


..C5J 





a-- 


XI IX XX 

12 12 22 

12 

A^^ = dot<AfB.O) = 


•«+ "a*-®*. 

®29 


®ax^ Sx'Sx 

s**- 

®;i2 


On expanding this we get 
^el = ^ol * Si •‘’22 **At .^2 


where : 


C63 


n: 


A = system open i™p clmrecUirisl-lf pilynomin] r ileUA) 

OL 

xx 


'6i 


= numerator polynomiel of tlw Xi--tn-6i transfer function 
obtained by applying Crnmer's Ruin to (2) 


~ numerator polynomiel of the Xit-to-Aa transfer function 

02 


N 


Xl- X2 

dl Sz 


H 


obtained by applying Cramer's Rule to (2) 
is defined as the coipling numerator and is given by 


XX X2 

'Ax Sz 


Sx Sa *ia 

Sx Sa **28 

Si Sa "33 
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From C4) : 


m ” V 

ti 11 


0.. 

a +b .G 

11 

13 12 22 

G_ 

a +b .G 

11 

22 22 22 

G.. 

a +b .G 

11 

32 32 22 


13 


33 


33 


Xi 

X3 

X3 



r b 

.G 

b 

.G 


11 

11 

12 

22 


b 

.G 

b 

•G 


21 

11 

22 

22 



.G 

b 

.G 


31 

11 

32 

22 



.C7) 


Solving for Xi/Xioi the transfer function function of Interestt 






a^ +b .G 

12 12 22 

*13 





•’« “ii 

a +b .G 

22 22 22 

*23 


Xi 

as « 


b .G 

31 11 

a +b iG 

32 32 22 

*33 


Xlo 


a ^b tG 

11 11 11 

a +b •G 

12 12 22 

*13 


• 

a 'll} .G 
. 21 21 11 

a +b .G 

22 22 ‘ 22 

*23 

: 


a +b '.G 
.31 31 11 

a +b iG 

32 32 22 

*33 


G„.CN!I? + G„.n!! 


'ir'"«i 


1 xZ 

zz'"6t 6z 


) 




..(8) 


Dividing the nunerator and the denoninator by A , and setting 

O 1 

iCVa = xi, 


'St. ol 


X3 


‘6i 

$2 


we get, 
Xi 




Xlo 
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22 


< & / N*; 


X 1 


) 

xl x2 


X 1 


1 + G«-»62 ^ i ; / 


)) 

rg^ 


Note tt«t this can be put in the form <r/(l '*' (')■ For this 
unity feedback closed loop systen, blie correcpcndiTig open Joo|.> 
system simply becomes G, given by : 




Cl + 




CIOJ 


Fig. 3 shoHS how the open loop transfer function has been 

modified by the two factors (If 632 -**61 62^61^ 

£1/(1 + ^ result of closing only the fe-to-^a loop. 

The new Xi-to-6t transfer function, can be rewritten as : 

. (» 6 ‘ * ^ > CUJ 

It 80 turns out that the zeros and poles of corre^>ond 

to the characteristic roots of the two rather siiple unity 

feedback systems shown in Fig. 4 a and Fig. 4 b respectively. The 

loop transfer functions of these two loops are, 

« 

reapeotively : 
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X 'O 

m ji 

X *0 

ss 

as 
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zz 
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Fig.3 Equivalent Block Diagram of the system : Xi — >6t, 





Fig. 4 b 


Fig. 4 Equivalent Sub-systems of the system : Xz 


Cl:\' t 


I :• T.. RAKpUfl 


• *'»!% I 
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-- and - 


, „x I M a 


M 


Xt 

6x 


^22 -”61 A2 


The above obeervation Is of nniclal sifinifioance as the pnlcs 
as Nell as the zeros of can now be roAuJaled by I.Ipj 

root- locus design of these two “feedbaok snb-systensi'' The 
suitable choice of G^j^now enables us to proceed with tl» final 
step of synthesizing Obviously this c*OBc>ensttitor oust be 

chosen so as to ensure tlie required overall system performance. 


APPLICATION OF THE PROPOSED DESK’iN APPHLIACH 
Tlie variables used In the synthesis are as fulLiws : 



Xt 


r 


iSiib 

X = 

Xa 

. »• . 

- 

u / II 

o 

; <!» = 



Since tlie primary objective is to ciailruL flight path, it is 
chosen . to the 'outer' feedback loop ami tpeed (u) in ctuson 
to form the 'inner' loop. 


The main problem being the non-minimum phase zero -1/T . . , we 

ni 

first try to place the zeros of tlie outer fliglit path loop by 
designing the closure of the coupling Loop B, since the closed 
loop poles of this one are tlie vory opMi Irnjii * 2 ^ros Mf tlie outer 
loop. Also, since the numerator rif tlm l.•>•l■t■r•sp(«ldirlg npni locE' 
transfer function i.e., tlie ixmplliifl luimvrnUir. Iia» iqin 
non-minimum phase zero, naturally It is essential tlial. Uiis loop 
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have no open loop poles in the right iKiJf plane. In ol.her worcJF:. 
we need to select so that it ha.«« no imn-miniRum pliase aerus. 

The only ohoioe available is tlie nunemtor of the ?■ -to-dab 
transfer function, and this fixes tlie choice of tlie two feedbacks. 

DESIGN DATA 

The following design data is assaiuetJ for n generic delta -wiiigeil 
high perfoniance fighter flying at sea level at. M 0.2, r = -2 deg : 
/ V = H0.514B ; g.sln = -0.00h0 ; - -0.2670 

=-0.87 ; =-1.85 ; M^.U = 0.1660 

/ U =-0.1450 ; g.oos = 0.1441 ; =-0.0765 

Z^/U = -0.1453 ; U = 0.0001 ; = -3.0050 

**dab "" -0.0015 ; X^^ / U r -0.0466 ; X^^/U = -0.0506 

The equations of notion of the aircraft are : 

-Mu-« 

®“\i <"V i U 

0 : u "IT 

-Z^ g.Bin g) (g.sinj; -s) o 

U U 

Here the state vector coeprisos of tho diuensionless forward 
velocity, the angle of attack and the flight path angle (instead 
of ttie pitch attitude angle normally used). 
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A|]|>lylng Cramer's Rule to the equations of notion we have 
the followii^ : 

1) Ihe open loop oharaoteristio polynomial i.e.i the denominator of 
the transfer functions is 

A a 8^ + 2. 4411. s^ ••• 1.9630.s^ + 0.1292.S + 0.0486 

O I 

a [s*+ 2C0.88)(1.36)s + (1.36)^][s‘+ 2(0.11)(0.i6)s + (0. 16)^1 


2) u N ” 
_ s 2f 

6* A 

ol 


- 0.0466 s^ - 0.0889 s^+ 0.8403 s + 0.2104 


^Ol 


- -0.0046 <s - 3.5424). Cs * 0.2449). (s + 5.2053) 

A 


ol 


3) T hJ. 


6m A 


ol 


0.1453 s^ + 0.2675 s^ - 1.4328 s 0.1200 


ol 

0.1453 (8 - 0.0852). (8 - 2.2964). Cs + 4.2226) 


ol 




5ab A 


ol 


0.0506 8^ - 0.1197 8^ - 0.0918 s - 0.0001 


ol 

0.0506 (8 + 0.0011). [s^+ 2(0. 88)(1. 35)8 + (1.35)^] 


^1 


5 ) 7 
5«b 


N I - 0.0001 ^ - 0.0136 8^ - 0.0258 8 - 0.0168 

Oab „ 


“ol ol 

-0.0001(8 + 135.12). l8^+ 2(0. 86) (1.12)8 * (1.12)^] 

A 


ol 
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6) Tbe coupling nunerator la given by : 

. 0.0074 8^ + 0.0136 a ' 0.0725 
a 0.0074 (a + 4.1824).C8 - 2.3446) 


N 7 « 

" 5«b 6m 


DESIGN 

Since tbe placenent of the outer-loop zeroa la of paramount 
luportanoe, aa alao explained earlier, we begin by deaignlng the 
oonpling loop 'flrat. The root looua technique la utilized here 
for the deaign. 

For the data aaauned. the qpen loqp tranafer function of the 
coupling 10 ( 4 ) can be written aa 

! 0.0074 a* + 0.0136 a - 0.0725 

^ 5«b - 0.0001 8^ - 0.0136 8^ - 0.0256 a - 0.0168 

Fig. 5 alxjwa the root looua plot for thia tranafer function. 
Aa can be a ee n . the preaenoe of the non~ainluuni phaae zero at +2.3 
tenda to draw the root looua Into the right half plane even 
for low valuea of gain. Thia will clearly not auit our purpoae. 

we try to remedy this by aeleoting a oonpensator that will 
pull the polea to the left whan the loop is closed. This required 
several trials, however. It was found that using the lead 



Ca ■+ 1.5) 

TT^nsr 


aooompl iahea the dealred task aultably. 
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Fig. 8 shONB the root locus plot with Uii?: i:OR(.>criRntiir in the 
loop. Selecting a gain of 0.3 we detain l.lw follfiwjng oJoenf1 loop 
poles : 

-0.8082 ± 0.3378 i ; -.•>.8364 ; -157.85. 

It nay be noted that since none of these is in tlie riglit half 
plane « our outer loop transfer function will have no non-niniMim 
phase zeros. 

With the coapensator now cltosen, tlie inner loop A is nlosed 
For this loop the open loop transfer function can be stated as : 

H ^ - 0.0468 s” - 0.0889 s% 0.8403 s f 0.2104 

s* + 2.4411. s” + 1.9839.S* 4 0.1292.S i- 0.0488 

The root locus plot for this tran^d'er function is sliown in 
Fig. 7 and the clcsed loop polos for the gain k = 0.30 are : 

-6.0033; -1.2217 ± 0.83881 ; -0.0042 ± 0.1313i. 

These are the open loop poles of the onu^r Imip. 

Having placed the open loop poles and zeros for the flight 
path loop, we now focus on designing the closure of tin outer 
loop. This is done using the 'connect' fimct.ion of the HATLAB 
software package which allows us to handle the 
lulti-input-nilti-output system at hand. The open loop transfer 
function of the outer loop obtained is : 

-0.0001s* - 0.0ie5s*- 0.1 155s*- 0.1563s - 0-®^. .. 

s®4 8.455 s*4 18.66 s®4 IJ.70 s*4 0.337. .s 4 0.197 



ng.7 Root Locus of Inner Loop A 

(wHh l93d rompen^ator) 
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= -g- groua ^ 15 7.85XR « ^ 

(s f e.003)[s* + 2<0.89>(J.3R^s * n.nn)*] 

[s^ + 2(0.92X0. §0 )iL+_C«L§®jll_ 

[8* + 2(0.03X0.13)8 + (0.13)*! 

and this clearly checks with the poles arnl eusifliied. 


Fig. 8 shows the root loriiB plot I'lii' l|ir> 'iveralL HIH() 
feedback control system resulting fruai tlie closure of this loop, 
fie select the gain k = 9 to ensure adetjuata ■ pltjgoid damping. 

The response of this closed loop system to a unit step }' 
connand is plotted in Fig. 9. He also obtain the response of 
the other important fli^t variable, airspeed, to tlie same input. 
(Fig. 10). 

.It may be observed that while the nirsiieed fluotuetiuris 
eventually vanish as required, tlioy do so after going tlirough a 
number of slow oscillations. It is therefore essential to go 
through anotlier iterative atteept -at designing tlie outer loop 
closure with a view to achieving guiok decay of the airspeed 
fluctuations. 

It is obvious that the lung decay time exhibited by tlie 
airspeed response -is a manifestation of the plugoid dynamics with 
its characteristic low frequency oscalluticti. Tlie cutest way U» 
reduce the decay tine would Ije to use a compensator that nullifies 
ttie effect of tlie plugoid poles. Hence after a fow trials, a second 



ng.8 Root locus of Flight Path Loc 





34 - 


order compensator is selected Hj1.h 

0 = (s^ 4 0 .tWO4 s 4- 

s * + S.Ti fs 4 3 

To check the improvenieni in tho systen perfomarice, the tine 

responses are again studied for unit strp inputs to (r/ 2 '^) 8 nd 

Cu/}' )• These are plotted in Fig.il and Fig. 12 resE'ectively. Ttie 
c 

remarkable effect of the compensator on outer loop response is 


iniediately apparent. 





.junotlJnviad 'nN iPW 


- 0.019 

“^.02 

- 0.021 
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CIIAP'l'ER 4 

RESULTS AND CONCLUSIUHS 

^ As oan be seen from Fig. Jl illustrating tlw controller 
offc^itiveness, a positive step Input Ua’ough the airbrakes results 
in a negative steady state flight path angle. This implies that 
extending the airbrakes nos results in o steady rate of descent. 
In addition there is no hint of a non ■ mlnlaun plisse beliavioni^ond 
the plot shows a well daaficd orir«p resprtise. 

The airspeed response for t)ie sane step infAit is presented in 
Pig. 12. Evidently there is a “svmll" dnip in tl» forwaid 
airspeed before it is quickly neutralised as tlie airplane desceiKls 
and accelerates. However when tlie steady state is atbained, there 
is no subsequent ohange in the airspeed. It may be emphasised that 
tiie resulting chmge in ttie aicspfwxl oLmorved is fairly small. 
Hem, the change in tl» flight path angle by much as 1 deg. 
results in a change in the airspeed by l«a: l.lmii 0.02 H w)iii*-h 
happens to bo insignificant even at Llm Inw qaeeds of the 
approecli. 

Thus, Uie desired objective of rectifying tlie anomalcms 
flight path response at low speeds has been aohlrved. E^tliermore, 
tiia control synthesis enables decoupling the fligiit path from tho 
airspeed. The importance of this cannot t«e overemphasised. It is 
now KiQSSible to select tho airopewl befor*? engaging the controller 
for approach. In view of tte deconplhig 8..‘hievod and with the 

<Jr«ilre^l 



38 


worry about the airspeed fluctuations. Ihus, the pilot is free to 
ronlrnl hfs flight path with virtually no concern regarding any 
iiiKlur rlwnges in tlm trimmed airspeed. Tills is clearly a much 
simpler procedure than the establisiied one of oontinuously 
JuSi!Hnf{ throttle and the stick. 

CX3NCIJUSI0NS 

A novel ai^roech at fliglit path control lias been attempted to 
farilHnle airplane handling during the landitig a|)proach. Jnftial 
emiilvisin was laid on theoretically establishing as to why the zero 
in tlie fliglit path - to -elevator transfer function is 
considerctJ as tlie parameter of iwime inqportance. Ihe oonceptml 
clarity achiex'ecl in tlie process proved to bo of considerable 
im|*irlamT in tlie <*onlroller syntliesls. It was realized tliat the 
plmM<mfnt of this zero lias to be regulated to ensure its location 
in tlM» left Imir s-plane with an overall view of achieving the 
dentred flight |«lli control. With tlw simple i-oot locus design 
iiHwept . niifdlc'l here to a mull 1 - Iniwit - imdli - output systimi, 
it iKis Iwcn iKjssihlc to achieve this browl objectlw comnendably. 


HRfWirdflWITONS FOR I'UlinHiR V/QHK 
lltr irfiwnt l» slwHaHc in imlure u It W>11« to n 

InnrtiK H)«d. It need, to be extended for geneml 
.IVtu-nlinm ewerliX »-er a r«n*c of Meoh nuabere- 

Me- rfferle of Kuet* ta'o nol bc-n ..•cmei.lrred. HeeK euel 

le- «ref..ll»- •» «»«l« '•'’o 

.kiruel Me- iwellim ewe«Kdi e* wll ee lie- flare - and - louehdoen 
inie.n>-r.HieiM>. Ibo .-ffw-le of ineludinS lh« airbrakes 



ni.'l.iinUir Ruinl. also be examitic^d. 


Oie nnn envisage a cnRir<W>1.e fliglih (.•nLh fiOiiLrullei- 


Uni: 


fcntally decouples flight path atvJ alrsi^cud fi-on cnoh utlier and 
includes all Jmxllng plwsos iiiuluiUiig flare aixl touclriowii . 
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